Surface Bound States in n-band Systems with Quasiclassical Approach 
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We discuss the tunnehng spectroscopy at a surface in muhi-band systems such as Fe-based su- 
perconductors with the use of the quasiclassical approach. We extend the single-band method by 
Matsumoto and Shiba [J. Phys. Soc. Jpn. 64, 1703 (1995)] into n-band systems (n > 2). We show 
that the appearance condition of the zero-bias conductance peak does not depend on details of the 
pair-potential anisotropy, but it depends on details of the normal state properties in the case of 
fully-gapped superconductors. The surface density of states in a two-band superconductor is pre- 
sented as a simplest application. The quasiclassical approach enables us to calculate readily the 
surface-angular dependence of the tunneling spectroscopy. 

PACS numbers: 74.20.Rp, 74.25.0p, 74.25.Bt 



I. INTRODUCTION 

Much attention has been focused on novel Fe-based su- 
perconductors since the recent discovery of superconduc- 
tivity at the high temperature 26K in LaFeAsOi_j;Fa;fi. 
Many theoretical and experimental studies on Fe-based 
superconductors have been reported for the last year. It 
is important to identify the superconducting order pa- 
rameter to elucidate the mechanism of superconductivity 
in those high-Tc materials. 

A ±s-wave pairing symmetry has been theoretically 
proposed as one of the candidates for the pairing symme- 
try in Fe-pnictide s up er conductors ^ ° i ^ ^ The 
lbs-wave symmetry means that the symmetry of pair po- 
tentials on each Fermi surface is s-wave and the relative 
phase between them is tt. Recently, we showed that a 
fully-gapped anisotropic ±s-wave superconductivity con- 
sistently explained experimental observations such as nu- 
clear magnetic relaxation rate and superfluid density. 

A key point to identify the ±s-wave symmetry is a 
detection of the sign change in the order parameters be- 
tween Fermi surfaces. It is difficult to detect the rel- 
ative phase of the order parameters in a bulk mate- 
rial. However, as shown in studies of high-Tc cuprates, 
Andreev bound states are formed at a surface or a 
junction when the quasiparticles feel different signs of 
the order parameter before and after scattering . ^'^'^"^i"'^^ 
Since one can extract the information on the relative 
phase through Andreev bound states, several theoreti- 
cal studies on junctions and surfaces have been reported 
recently.^>ii'i^i2ai2ii2^i2^i2i Andreev bound states at 
zero energy have been experimentally observed as a zero 
bias conductance peak (ZBCP) in tunneling spectroscopy 
for Fe-based superconductorsi^ 

The Fe-based superconductors are interesting also as 
novel unconventional multi-band superconductors since 
multi-band effects are essentially important therei^ Fermi 
surfaces in these systems predominantly consist of the 



d orbitals of Fe atom. Kuroki et al.^ suggest that five 
orbitals are necessary to describe the properties of the 
superconductivity, and they elaborate an effective 5-band 
model. On the other hand, MgB2 is a 2-band system that 
is a conventional s-wave BCS-like superconductor. 



The aim of this paper is to develop a method for ana- 
lyzing surface bound states in multi-band superconduc- 
tors. Matsumoto and Shiba^^ developed a method to 
analyze surface bound states in single-band systems such 
as high-Tc cuprates. We extend their method into multi- 
band systems. Since the ratio of the superconducting 
gap A to the Fermi energy Ep is small, A/Ep <C 1, 
in Fe-based superconductors, we can adopt a quasiclas- 
sical approach. In this approach, all we need is only to 
consider quasiparticles at the Fermi level. Thus, we can 
reduce computational machine-time and the physical pic- 
ture becomes clear. In addition, this approach enables us 
to easily calculate the surface-angular dependence of tun- 
neling spectroscopy. We find a general appearance con- 
dition of the ZBCP for multi-band systems. This general 
condition can be applied to various pairing symmetries 
including its-wave and d-wave. With our method, we will 
discuss a two-band superconductor as a simple example. 



This paper is organized as follows. The formulation of 
our quasiclassical approach is shown in Sec. II. We ap- 
ply a quasiclassical approximation to eliminate fast spa- 
tial oscillations with Fermi wave length. The appearance 
condition of the ZBCP in multi-band systems is derived 
in Sec. III. The results for a two-band model are shown 
as a simple example of our approach in Sec. IV, where 
we will show both analytical and numerical results. The 
discussions and conclusion are given in Sees. V and VI, 
respectively. In the appendix, we describe the derivation 
of the appearance condition of the ZBCP when a system 
can be treated without quasiclassical approximation. 
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FIG. 1: Schematic figure of a specular surface. 

II. FORMULATION 

A. Orbital representation and Band representation 

Let us consider the local density of states near a sur- 
face following a procedure by Matsumoto and Shibai^ 
We assume a two-dimensional superconductor, and con- 
sider a specular surface, for which the component of the 
quasiparticlc momentum along the surface is conserved 
as shown in Fig. [T] We treat the surface as a potential 
t/(r)f3, where the time-reversal symmetry is conservedi^ 
Here, fi [i = 1,2,3) denote Pauli matrices in Nambu 
space and r is the position in the real space. We con- 
sider a rt-orbital system, which is a periodic crystal with 
n atomic orbitals in unit cell. Throughout the paper, hat 
a denotes a, n x n matrix in the orbital space, and check 
a denotes a 2n x 2n matrix composed of the 2x2 Nambu 
space and the nxn orbital space. We calculate the Green 
fimction under the influence of U{r)f^. It is written as 

G(r,r') = Go(r,r') + j dr" Go{r ,r")U {r")f:,G{r" ,r). 

(1) 

Here, Gq is an unperturbed Green function in the absence 
of U . We take the a;(y)-axis perpendicular (parallel) to 
the surface as shown in Fig. [T] Considering the surface 
situated at x = and the scattering potential U written 
as U{r) = Uq6{x), Eq. ([1]) is reduced to 

G{X, ky,x', k'y) = 27rSiky - k'y)G{x, X' , ky), (2) 

where 

G{x,x',ky) = Gq{x,x' ,ky) + Goix,0,ky)Uof3 

X (1 - Go(0, 0, fcy)C/of3)-^Go(0, x', ky). 

(3) 

Here, we have taken the Fourier transformation with re- 
spect to y. We use units in which h — 1, and the coor- 
dinates r and the momentum k are dimensionless. The 
surface is actually represented in the limit Uq oo. The 
Green function is then given by 

G{x,x',ky) = Go{x,x',ky) + Gp{x, x' , ky), (4) 

where 

GP(X, X', ky) = -G0{X, 0, ky)Go{0, 0, ky)-^Go{0, X\ ky). 

(5) 



The local density of states at the position x for the mo- 
mentum ky is written as 



N{x,ky)^ Im [TrG^(a;,a;,fc„)], (6) 

TT 

where 



G [X, X, ky) G[X, X, fcy)|i(^^_,£;_|_j^. 



(7) 



Here LUm is the fermion Matsubara frequency and 77 is a 
positive infinitesimal quantity. The unperturbed Green 
function G^{x, x' ,ky) is given by 

G^ix,x',ky)^^ f dk,e'''^^^-^'^G^{k,,ky), (8) 



where 



(kx, ky) — [E — H^{kx, ky)) 



(9) 



Here, H^{kx,ky) is the 2n x 2n Hamiltonian in Nambu 
and orbital spaces written as 



H° A° 



(10) 



in the "orbital representation" where the base functions 
are atomic orbitals in crystal unit cell. From now on, 
the subscript "o" indicates that matrices are represented 
with the orbital basis. H° is the Hamiltonian in the 
normal state represented as n x n matrix in the orbital 
space. Remember that n is the number of the orbitals. 
A° is the superconducting order parameter. 

Let us introduce a n x n Hamiltonian in the "band 
representation" defined by 

H (kx^ky) — P {kx,ky)H [kx, ky)P[kx, ky), (11) 
Ai 

■■. I- (12) 
A„ 

Here, Ai (i = 1, 2, • • • , n) denote the eigenvalues where 
the relation A; > A^ {i < j) is satisfied. P is a unitary 
matrix consist of the eigenvectors that diagonalizes the 
Hamiltonian H°. The 2n x 2n Hamiltonian in Nambu 
and orbital spaces in the "band representation" is also 
defined by 

H^{kx,ky) — U [kx, ky)H-j^[kx, ky)U{kx, ky), (13) 



(14) 



where 



A'' = p-iA°P. (16) 
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In general, A'' contains ofF-diagonal elements, which cor- 
respond to inter-band pairings. Assuming that intra- 
band pairings are dominant, we neglect the off-diagonal 
(inter-band) elements in A**: 



/ Ai 

■■. 

V A„ 



(17) 



That is, we consider that only single pair-potential is 
defined on each Fermi surface. Here, A^ is the pair- 
potential on the i-th band. Substituting Eq. (|13p into 
Eq. the Green function G^{kx, ky) is written as 



(18) 



Assuming Eq. (|17p and taking the inverse matrix of — 
H'', one can obtain 



A+ B 



(19) 



where 



B = 



-|Ai|2+£;2_Aj " 

V 

/ 
v 






Ai 

-|AiP+£;2-A? 














A„ 



^ -|A„|2+£2-A2 / 

We find that Eq. can be rewritten as 

G^ikxi^y) — ^ ^ G* (fca; fey); 

i 

where i is the band index and 
1 



(20) 



(21) 



(22) 



[M, 



(£; + AOM, A,M, 
-|A,;|2 + £2 - Af V A*M, (^ - A,)M. 




2nd channel (/=2) 
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E„ 
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FIG. 2: Schematic figures of band dispersions along a fc^, line 
witli a fixed k^. 



and Fe-based ones. In this case, one can use a quasiclas- 
sical approach. 

We consider a line with a fixed ky in the momentum 
space. On this line, we classify n-bands into two groups. 
One group is composed of the bands on which the eigen 
energy Xi^k^, ky) crosses the Fermi level (for example, the 
bands i = 1 and 2 in Fig. [2]). The other group is com- 
posed of the bands on which the eigen energy does not 
cross the Fermi level (the band i = 3). For the former 
group, we can analytically integrate &(kx,ky) over kx 
with the use of a quasiclassical approach since G^{kx, ky) 
is a function localized near the Fermi level. For the latter 
group, we need to integrate G^{kx,ky) over kx numeri- 
cally since G^{kx,ky) is not a localized function. How- 
ever, the integrand is a smooth function, so that it is easy 
to perform such a numerical integration. 

We integrate &{kx, ky) on the bands of the first group 
with the use of the quasiclassical approach. To perform 
the fcjj-integration, we divide the kx-Vme with a fixed ky 
into some segments as shown in Fig. [51 Each segment 
has only single channel that is the point satisfying the 
relation Ai = Ep. From now on, I denotes the channel 
index and k denotes the maximum number of I. The 
integration for the i-th band is written as 



j-[J-ooV^{X^) 



(26) 



(23) 
(24) 



Equation ([5^ is divided into a sum of the Green func- 
tions defined on each band. Substituting Eq. ((22)) into 
Eq. ([8]), Gq{x,x' ,ky) is expressed as 

G^{x,x',ky)^Y.^ J dkxe'''''^^-^'^G\kx,ky). (25) 

i 

Hence, the fc2:-integration is found to be performed on 
each band independently. 



B. Quasiclassical Approach 

We assume |Ai| ^ Ep. This relation is satisfied in 
most of systems such as conventional superconductors 



Expanding kx{Xi) in the first order of Ai around Ai = 
Ep as kx{Xi) = kpx + Xi/vpx, one can carry out the 
integration by the residue theorem: 

^ j dkxe^''^^^-^'^G\kx.ky) = -^j2GfMx), (27) 

1=1 

where 

^fA^Fx) 



f±{kl 



f+ (fc>L)M. (fc>L) A. (fck)M. {k'/J 

A,(fc>'jM,(fc>'j f.{el)Mel) 



Hk^x). (28) 
I 



(29) 



E ± sgn (x - a;')sgn («>L)V^' " 



(30) 
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Here, kp^ and Vp^ are the Fermi wave number and the 
Fermi velocity on the Z-th channel, respectively. Using 
the above, Eq. ([25|l can be written as 



Gf){x, X ,ky) — 



ieQ 1=1 



Fx) 



+ T.i j dk.,e'^^^--^'^G\k,,,ky), 

(31) 

where the elements in Q are the indices of the bands 
whose energy dispersions cross the Fermi level for a fixed 
ky. Here, we assume Aj^g — 0, namely the super- 
conducting order parameters are finite only around the 
Fermi level. It should be noted that the second term 
in the right-hand side of Eq. (|5T|) cannot be neglected 
since G^{^,Q^ky)~^ without this second term may have 
artificial divergences. 



C. Eliminating the fast oscillations with Fermi 
wave length 

We assume the condition kp^ 3> 1 (i.e., |Ai| ^ Ep), 
which is the quasiclassical condition. Here, is the coher- 
ence length of a superconductor. Under this condition, 
the short range spatial oscillations characterized by the 
Fermi wave length Xjkp can be eliminated. We rewrite 
Eq. (01]) as 



G'g '(x, X , /Cy) 

where 
KieQ{kj:,ky) = 

Ki<^Q{tx,ky) = 



^ ^ J dkxKi[kx,ky)e 



ikx {x — x ) 



(32) 



\Y,Gl{x,x\kx)5{kx-ej}^), (33) 



—G^{kx, ky). 
Ztt 



(34) 



The perturbed Green function Gp{x,x',ky) defined in 
Eq. (IH) can be written as 



G^{x , X , ky) 



J2 J dkxdky^''='^-''>">K,{kx,ky) 



xG^{0,0,ky)-^k,„ik':,ky). (35) 

Setting exp[i{kxX — k'^x')] 1, we eliminate the short 
range oscillation while keeping the enveloping profile of 
the integrand. Thus, the above equation is reduced to 



Gf p{^x ^ X , ky) 



J2 J dkxK,{kx,ky)G^iO,0,ky)~' 

i 

J2jdk'^K^"{k'^,ky). (36) 



From this equation, it is concluded that the Andreev 
bound states appear when Gq {0,0, ky)^^ diverges, i.e., 
when det G^ {0,0, ky) = 0. 



III. APPEARANCE CONDITION OF THE 
ZERO BIAS CONDUCTANCE PEAK (ZBCP) 

Let us consider the appearance condition of the ZBCP 
in n-band system at a surface. At the zero energy E ~ 0, 



G(i{x = 0,x' = 0,A:^;) defined in Eq. 
written as 



2B [for ieQ] is 



sgn (Ai 



2\v 





M,{k'-^ 



Fx) 



M^{k'^x) 




For i ^ Q, we have from Eq. ((23| with E ^0, 



G' 



-A 



A,Af, 
-\,M. 



(37) 



(38) 



where we have set A; — because the superconducting 
order parameter is assumed to be finite only near the 
Fermi level and the bands with the indices i ^ Q do not 
cross it. Substituting the above equations into Eq. (|3ip . 
we can obtain the appearance condition of the ZBCP 



from det G^{0,0,ky 



det 



0: 

-/ L 

L i 



= 0, 



where 



leQ I "^l^Fxl 



I = 



E- 



27r y Xi{kj:) 



M,{kx). 



(39) 

(40) 
(41) 



Equation (|40p shows that the appearance condition does 
not depend on the anisotropy of the pair-potentials and it 
depends only on the signs of them because information on 
the pair potentials is included in the form, sgn {/S.i{k^p ^)) , 
in Eq. (j40| . This result shows that information on the 
normal state (i.e., the matrices Mi, v'p^) is important for 
the ZBCP to appear. 



IV. TWO-BAND MODEL AS A SIMPLE 
EXAMPLE 

A. Model 

We calculate the density of states in a two-band super- 
conductor as a simple example. We consider a two-band 
tight-binding model on a square lattice. There are two 
orbitals on each lattice site. The Hamiltonian with a 2 x 2 
matrix form in the normal state is described as 



^° ^ I 2t' sin 



—t cos(fca) — /i 2t' s\n{ka) sin(/cb) 
^' '''in(fca) sin(A:fc) —t cos{ki,) — /i 



(42) 



in the orbital representation (n = 2). Here, ka and kf, are 
the axes fixed to the crystal axes in the momentum space, 
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(a) 




(b) 




FIG. 3: (Color online) Fermi surfaces in the two-band model, 
(a) the half filling {fi = 0) and t' = O.lt. (b) = 0.2t and 
t' = 0.2t. 




FIG. 4: (Color online) Fermi surfaces — Q and t' = O.lt) 
and kx line with fixed ky. 



t and t' are intra- and inter-orbital hopping amplitudes, 
respectively, and ^ denotes the chemical potential. We 
use the unit in which the lattice constant a = 1. This 
Hamiltonian can be diagonalized into the matrix in the 
band representation, i?**, written as 









As 



(43) 



Here, Xa{b) denotes the energy dispersion on the A{B)- 
band. As shown in Fig. [3l the Fermi surfaces consist of 
two parts near the half filling. 

We consider the two-band s-wave superconductor de- 
scribed by the pair potential in the band representation: 





Ab 



(44) 



Here, A^(b) is the pair potential on the A{B)-ha,nd. 

We introduce the coordinates {ka,kb) fixed to the crys- 
tal axes: 



cost 
sin t 



— sin 9 
cos 9 



(45) 



Here, the kx{ky) axis is the axis parallel (perpendicular) 
to the surface and 9 is the angle between the fc^ and k^ 
axes. Considering [110] surface, we fix 61 = 7r/4. The 
quasiparticle momentum ky is conserved since we con- 
sider the specular surface. 

It should be noted that one needs to treat the Bril- 
louin zone in the surface-coordinates {kx,ky) for each 
surface angle since it is necessary to consider all possi- 
ble scattering processes at the specular surface (namely, 
all fcy-momentum conserving processes). For example, it 
naively seems iir Fig. 3(a) that possible scattering pro- 
cesses occur only on the inner Fermi surface (red) for 
the [110] surface {9 = n/4) in the region 7rV2/4 < 
ky < 7r-\/2/2 [the ky axis is directed in the direction 
of {ka,kb) = (-1,1) in Fig. 3(a)]. However, for [110] 
surface, one has to consider also the outside of the first 
Brillouin zone as shown in Fig. 4 so that scattering pro- 
cess between outer Fermi surface (green) and inner Fermi 
surface (red) can occur. 



At the half filling for [110] surface, the second term in 
Eq. (f3T|) does not exist since the energy dispersions of the 
A and B bands always cross the Fermi level on k^ line 
with any fixed ky in the momentum space as shown in 
Fig. m In this case, / defined in Eq. (|4ip is zero because 
there is no band with the index i ^ Q. Therefore, the 
appearance condition of the ZBCP in Eq. (|39p can be 
rewritten as 

det L = 0. (46) 
where L is defined in Eq. (|40p . 

B. Analytical Results 

1. At the half filling for [110] surface 

We will analytically show that the ZBCP always ap- 
pears for any strength of the inter-orbital hopping t' in 
the case of [110] surface {9 — 7r/4) at the half filling. 
On the lines which satisfy ka — {kx — fcj,)/\/2 = titt or 
kb = {kx + ky)/^/2 = mr in the momentum space, one 
can easily obtain the unitary matrix P that diagonalizes 



P{ky) = <^ 

















I) 



ikx ky — ^oTT, or kx ^ ky — ti^^tt, 

(47) 



Here, ne(o) is an even (odd) integer. Substituting these 
P{ky) into Eq. (gHl), we obtain L: 



L cx 



sgn(AA) +sgn(AB) f 1 
1 



\VFx 



(48) 



The appearance condition of the ZBCP [Eq. (US])] is writ- 
ten as 



sgn(A^) -)-sgn(AB) = 0, 



(49) 
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on the lines where ka = {kx — ky)/^/2 = nn or kh = 
{kx + ky)l\/2 = nn. This condition is always satis- 
fied in the sign-reversing s-wave (±s-wave) supercon- 
ductors in this model. The its-wave symmetry means 
that the symmetry of pair potentials on each Fermi sur- 
face is s-wave and the relative phase between them is 
^^^^^5^6,7,8,9,10^11^12 Therefore, the ZBCP appears at the 
points on the Fermi surfaces where the relation ka = 
or fcfc = is satisfied in the momentum space. 



2. Case of t' /t = for [110] surface 

In the case of t' /t — 0, we can analytically show that 
the ZBCP always occurs for [110] surface. In this case, 
the unitary matrix P can be written as 



(a) al Half filling {]i = 0) 



P{ky) = 



1 
1 

1 

1 



kxky > 0, 



ky "\ . 



(50) 



As in the case of Eq. ([47|) , these unitary matrices lead to 
the same appearance condition of the ZBCP as Eq. 




FIG. 5: (Color online) The density of states at the surface 
for various surface angles. The pair potentials are Aa = A = 
O.OOlt and Ab = -Aa. (a) the half filling (/i = 0) and (b) 
fi = 0.2t. The inter-orbital hopping amplitude is t' — O.lt. 
The smearing factor is ry = 0.1 A. 
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3. Case oft'/t = at the half filling for [110] surface 

Finally, we discuss the difference between the appear- 
ance conditions with and without the quasiclassical ap- 
proach. As shown in Appendix, the appearance condition 
obtained without the quasiclassical approach for t' = 
at the half filling is written as 



/i = or 



0, 



0, 



(51) 
(52) 



where 



h. 



In 



(sin(±fc„/V2) + Vl + |A^/tP)" 
(sin( ± fe„ / V2) - Vl + IAlTTj^) 2 



In 



(sin(±fc^/v/2)-H^l + |AB/tp)" 
(sin(±fc„/y2)-^l+|AB/i|2)2 



2v/1 + |AbA|2 
sgn(AB/i) , sgn(As/<) 



v/i + |AaAP v/i + |AbAI^ 



(53) 



.(54) 



Here, we assume that the pair-potentials A^ and A^ 
do not depend on k for simplicity. The above equations 
suggest that the appearance condition of the ZBCP de- 
pends on the details of the amplitudes |A/i| and |Ab| 
in contrast to the quasiclassical result [Eq. (|15)) ] . In the 
limit of |Ayi_B/i| <^ 1, on the other hand, Eqs. ((53)) and 
([54|l are reduced to Eq. (j49|l obtained by the quasiclas- 
sical approximation. Thus, the quasiclassical and non- 
quasiclassical results coincide in this limit. Therefore, it 



FIG. 6: (Color online) The density of states at the surface 
for various surface angles. The pair potentials are Aa = A = 
O.OOlf and As = -0.5Aa. (a) the half filling (/i = 0) and (b) 
fj, — 0.2t. The inter-orbital hopping amplitude is t' = O.lt. 
The smearing factor is rj = 0.1 A. 



is suggested that our quasiclassical approach is appropri- 
ate when \AA,B\/t 1. 



C. Numerical Results 

The density of states at the surface is calculated from 
Eq. (O as 

N{E) = 7^ / dkyN{x = 0,ky). (55) 
2tt J 

We consider the is-wave 

superconductori2i^!^i^i^ii2ii2iiiS'iiii^ and the same two- 
band model as discussed in this section. 



1. Dependence of the surface-angle 9 

We show the energy dependence of the density of states 
for various surface-angle 9 in Figs. and El The peak po- 
sitions of the Andreev bound states depend on the surface 
angle 9. By comparing the results between Figs. O and 
ini it is noticed that those positions do not depend on the 
pair-potential amplitude. 
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(a):at Half filling (n = 0) {b):at |i = 0.2t 



I 1 1 1 1 





i. t' = • 




tl = 0.1t 




tj = 0.2t 




^ : 
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FIG. 7: (Color online) The density of states at [110] surface 
for various inter-orbital hopping amplitude t' . The pair po- 
tentials are A a = A = O.OOlt and As = — Aa- (a) the 
half filling (/i = 0) and (b) = 0.2t. The smearing factor is 
rj = O.IA. 



2. Dependence of the inter-band hopping amplitude t' 

We investigate the dependence on the inter-orbital 
hopping amplitude t' . We consider [110] surface {9 = 
7r/4). As shown in Fig.[7{a), the ZBCP always exists at 
the half filling (/i = 0) for any inter-band hopping ampli- 
tudes t'. At /i = 0.2t as shown in Fig. [Tljb), the ZBCP 
only appears when without an inter-band hopping, i.e., 
t' = 0. These ZBCPs appear when the appearance con- 
dition in Eq. (j49| is satisfied. 



V. DISCUSSION 

The advantages of our method are that one can eas- 
ily investigate the surface-angle dependence of the den- 
sity of states with the use of the quasiclassical method 
and easily calculate the density of states in the n-band 
(multi-band) system with less computational machine- 
time. Therefore, we can take, for example, a realistic 5 
band model in order to discuss the density of states for 
iron-based superconductors. We will report its results 
elsewhere near future. 

We have assumed that the matrix of the pair poten- 
tial in the band-representation does not have off-diagonal 
elements, which correspond to the inter-band pairings. 
When the inter-band pairing is dominant, the Cooper 
pairs have center-of-mass momentum q ^ 0. Usually 
such pairs are not energetically favorable since the pair 
potentials have spacial dependence even in bulk systems. 

Starting with the same Matsumoto-Shiba method, 
Onari et alM- recently calculated the surface Andreev 
bound states without the quasiclassical approximation. 
Their results show that the peak positions of the Andreev 
bound states depend on the gap amplitudes on two bands 
in the same two-band model as considered in Sec. IV, 
and the ZBCP does not always appear at the half-filling. 
These results might seemingly be inconsistent with our 
quasiclassical results. It is, however, not the case. 

They obtained the perturbed Green function by di- 



rectly integrating the original unperturbed Green func- 
tion over kx and ky numerically. The original unper- 
turbed Green function has sharp peaks on Fermi surfaces 
in the momentum space and rapid Fermi-wave-length os- 
cillations in the real space. We have integrated out those 
properties by the quasiclassical approximation. It should 
be noted that the pair potentials are of the order A ^ O.li 
in Ref. 19. This parameter is out of our quasiclassical ap- 
proach (A/t ^ 1). As shown in Sec. IV. B. 3, our analyti- 
cal result, which depends on the details of the gap ampli- 
tudes and therefore is consistent with Ref. 19, is reduced 
to the quasiclassical result in the limit A/t ^ 1. Thus, 
the differences in the obtained results between Onari et 
alM and the present paper would be due to the difference 
in applicable parameter regions. 

The formulation derived in Sees. II and III can be ap- 
plied to general multi-band superconductors including d- 
wave pairing superconductor. The appearance condition 
for the ZBCP is given as Eq. ^ in Sec. III. Let us 
consider, for instance, the case of the two-band model 
discussed in Sec. IV at the half filling for [110] surface. 
From Eq. ([M)) . the appearance condition for the ZBCP 
is given as 

sgn(A^i)-f sgn(A^2) + sgn(ABi)-Ksgn(AB2) = 0. (56) 

Here, A^i and Aa2 are the pair potentials on the inner 
Fermi surface (red) in Fig. 21 and A^i and Ab2 are the 
pair potentials on the outer Fermi surface (green) there. 
For a two-band d-wave superconductor, A^n = —Aa2 
and Abi = —Ab2, so that the above condition is sat- 
isfied and the ZBCP appears. Furthermore, in the case 
of a single-band d-wave superconductor, the appearance 
condition for the ZBCP is obtained from Eq. ([55)1 as 

sgn(AAi) +sgn(AA2) = 0. (57) 

This is consistent with previous results for d-wave pairing 
in Refs. [Ullililii, where the ZBCP appears when the 
quasiparticles feel a superconducting phase change tt in 
the surface scattering process A1^A2 on a Fermi surface. 

VI. CONCLUSION 

In conclusion, we extended the single-band method 
by Matsumoto and Shiba^^ into general n-band case 
{n > 2). With the use of the quasiclassical approxima- 
tion, we developed the way to integrate the unperturbed 
Green function with respect to k^ which is the momen- 
tum component perpendicular to a surface. We showed 
that the appearance condition of the ZBCP does not de- 
pend on any anisotropy in the pair-potential amplitude, 
but only on the relative phase, in the case of A ^ Ep in 
rt-band systems. The properties of the normal state are 
influential for the ZBCP to appear. 

We also calculated the surface density of states in the 
two-band system as a simple example of our approach. 
We suggested that our quasiclassical approach is appro- 
priate when \A\/t ^ 1. We showed that the peaks of 
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the density of states due to the Andreev bound states 
depend on the surface angle and the parameters in the 
normal state {t,t', fi), so that the sign-reversing s-wave 
(±s-wave) superconductors exhibit complicate properties 
in the tunneling spectroscopy compared with single-band 
d-wave superconductors. 
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where 



A, 



B = 



E±ei 






E±e2 






To investigate the appearance condition of the ZBCP, we 
set E = (i.e., zero energy), a; = and x' = Q (i.e., at 
the surface). Then, we calculate Gq{E = 0,x = 0,x' = 
0,ky): 



G^{E = 0,a; = a;' = 0,ky) 



dk 



V)- 



(A.8) 



Each element in this matrix can be integrated analyti- 
cally as 



APPENDIX: INTEGRATION WITHOUT 
QUASICLASSICAL APPROXIMATION 

We will show the ZBCP appearance condition 
Eqs. (jSip and (j52[) for zero inter-orbital hopping am- 
plitude {t' — 0), by integrating Eq. (P5|) in the simple 
two-band model for [110] surface at the half filling^Z 

The Hamiltonian in the normal state is described as 



ei{kx,ky) ^ 
e2{k^,ky) 



with 



^l{kx:ky^ — t COs{kx ^y): 

f-2{K,ky) = -tcos{kx + ky). 



(A.l) 



(A.2) 
(A.3) 



Here, we have introduced kx = kx/\/2 and ky — ky/^/2. 
Considering the pair potentials Aa,b which do not de- 
pend on k and using the unitary matrix Eq. (jSOp . the 
pair potential matrix in the orbital representation can 
be written as 



A" 



AAe{kx) + Aeei-kx) 

AAe{k) + ABe{-kx) 







(A.4) 



The unperturbed retarded Green function Gq {E, k^, ky) 
is written as 

G^iE,kxrky) = {E-H?,r' = / )(A.5) 



\A\^ + cos^{kx±ky) 
\A\^ + cos^{kx±ky) 



—dx 



dx 



|A|2 



(A.9) 
(A.IO) 



Integrating Gq{E = 0,kx,ky) by using the above formu- 
lae, we finally obtain Gq{E = 0, a; = 0, x' = 0, ky): 



G^{E = 0,x = 0,a;' = 0,ky) cx 



-h 





Aab 













Aab 


Aah 





h 








Aab 





I2 



(A.ll) 



where /i_2 and Aab are defined in Eqs. ([55]) and ([51]) . Its 
inverse matrix is written as 



/ 



cx t 



A„6 n h 



\ 















(A.12) 



The zero energy bound states appear when [G^]"^ di- 
verges as noticed from Eqs. dU and ([S]). Therefore, the 
appearance condition of the ZBCP is expressed as 



Aab - 0, 
/i = or /2 = 0. 



(A.13) 
(A.14) 
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